Informational completeness of continuous- variable measurements 
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We justify that homodyne tomography turns out to be informationally complete when the number of indepen- 
dent quadrature measurements is equal to the dimension of the density matrix in the Fock representation. Using 
this as our thread, we examine the completeness of other schemes, when continuous-variable observations are 
truncated to discrete finite-dimensional subspaces. 
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Introduction. — Measurement lies at the very heart of quan- 
tum information. It is an indispensable tool for identifying 
how well one can prepare or create a particular quantum state, 
as such an assessment is achieved by performing suitable mea- 
surements on a sequence of identically prepared systems. A 
set of measurements whose outcome probabilities are suffi- 
cient to determine an arbitrary quantum state is called in- 
formationally complete (IC) 111, while the process of recon- 
structing the state itself is broadly called quantum tomogra- 
phy d. 

Investigations on IC measurements have been extensively 
carried out for discrete Hilbert spaces [3|]. For continuous 
variables, the archetypical example is homodyne tomography, 
first suggested in the seminal paper of Vogeljind Risken ||4[] 
and implemented by Raymer and coworkers JH]. Still, the re- 
construction of the measured state is inexorably accomplished 
in some finite subspace chosen by an educated guess. 

In many experimental situations one has some additional 
prior information that can be efficiently used for reducing the 
number of relevant unknown parameters, making the recon- 
struction much simpler. For example, the knowledge about 
the Gaussian character of the state reduces the problem to the 
evaluation of its 2 x 2 covariance matrix, which is formally 
analogous to the estimation of a spin 1/2 In the same vein, 
if the system can be represented (say, in the Fock basis) by a 
finite density matrix, the number of quadratures needed for an 
accurate reconstruction in homodyne tomography is precisely 
equal to the dimension of that matrix 01- Note, in passing, 
that this is indeed very closely related to the quantum version 
of the sampling theorem [HI]. 

With the advent of powerful nonlinear estimation tech- 
niques, offering amazing performance and robustness in most 
applications, matching the signal space used for coding the in- 
formation in an experimental setup is more relevant than ever. 
This is, in fact, a dual problem: one is not only interested in 
what can be reconstructed from a particular scheme, but also 
in identifying the best option for a given signal. Any mea- 
surement is represented by an operator and what matters for 
tomography is the number of linearly independent operators 
required for the reconstruction. Keeping in mind the lesson 
from homodyne detection, the question of how to select those 
operators appears far from trivial. 

In this work, we examine the informational completeness of 



various schemes. We derive a surprisingly simple connection 
between the different measurement settings and the number 
of linearly independent elements induced by them. In addi- 
tion, the analysis displays the information provided by each 
acquired piece. Such knowledge is certainly of great interest 
both for the optimization of feasible tomographic protocols 
and for the design of new ones. 

Preliminaries. — A cZ-dimensional quantum system is rep- 
resented by a positive semidefinite d x d density matrix that 
requires d 2 — 1 independent real numbers for its specification. 
A von Neumann measurement fixes at most d — 1 real param- 
eters, so d + 1 tests have to be performed to reconstruct the 
state. This means that d 2 + d histograms have to be recorded: 
the approach is thus suboptimal for this number is higher than 
the number of parameters in the density matrix. The von Neu- 
mann strategy can be further optimized regarding this redun- 
dancy when the bases in which the measurements are per- 
formed are mutually unbiased H]. 

Besides, it is known that more general measurements ex- 
ist. Such generalized measurements appeared as positive 
operator-valued measures (POVMs) in the quantum theory of 
detection in the early 70s. It was soon realized that, for nu- 
merous reasons, they provide greater efficiency flol - Hill . 

In short, a POVM is described by a set of linear opera- 
tors {E„} furnishing the correct probabilities in any experi- 
ment (we assume for simplicity discrete outcomes) through 
the Born rule 



p n = Ti(pE n ), 



(1) 



for any state described by the density operator p. Here Tr 
denotes the trace of a complex matrix. Compatibility with the 
properties of ordinary probability imposes the requirements of 
positivity, Hermiticity and resolution of the identity: 



E n >o, e„ = eI ££„ = i. 



(2) 



The case of projective von Neumann test is recovered when 
all the operators in the set {E n } commute. 

A quite general way of implementing a POVM is to prop- 
erly combine the system to be tested with an auxiliary system 
whose initial state is known and then to perform a von Neu- 
mann test on the combined system lfl5ll . When the coupling 
with the auxiliar system and the von Neumann measurement 
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are judiciously chosen, one can get IC POVMs, which in ad- 
dition are optimal in the sense that minimize the number of 
independent detections that must be obtained during the tomo- 
graphic process |16fl. Below, we shall exploit a similar idea to 
make IC measurements from incomplete ones by cutting the 
dimensionality of the signal space IU7I1 . 

Quadrature measurement. — To get the backbone of our 
proposal, we start with states that can be written as a finite sum 
in the Fock basis. We wish to characterize a ^-dimensional 
subspace i ( / = ]Cn=o \ n ) ( n \ from projections onto the quadra- 
ture eigenstates \xg), which can be expressed as 

(n\x e )ocH n (x)e-^/ 2 e ine . (3) 

Our task is find out how many different phase settings (also 
referred to as cuts) are needed to make the scheme complete. 
It is clear that the complex argument of (n|jce) is indepen- 
dent of x, so one quadrature can never generate a complete 
POVM. For example, the three states (|0)(0| + |l)(l|)/2 and 
(|0) ± i\ l))/v2 have the same position distribution and hence 
the position measurement (0=0) cannot discriminate among 
them. Obviously, since (n\x$-o) is rea l f° r a U *> sucn projec- 
tions fail to generate the oy component. This difficulty persists 
in higher dimensions. 

Fortunately enough, this topic can be fully analyzed in a 
closed form. Indeed, the probability distribution of the homo- 
dyne detection can be written, up to a normalization factor, 
as 

p(x, 9) « £ Pkl H k {x)H t {x)e-* V^ e , (4) 
k,e=o 

where x is the measured amplitude. Let us take all matrix 
elements non-zero and calculate the number N of linearly 
independent combinations of the density matrix © that can be 
generated from different choices of a: and 6. For every poly- 
nomial term with a maximum power x", we have the following 
different combinations of Hermite polynomials [omitting for 
brevity the associated factor e l ( k ~ ' 6 ~\: 
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(5) 



For a single quadrature 0\ , the number N\ is determined 
by the term Hd-\Hd-\ with highest polynomial power and 
is equal to A^i = 2d — 1 (the total number of first elements 
H\fi\ in all rows of the above array). Since we have all the 
polynomial powers from up to 2d — 2, there are no more 
linearly independent elements. The second value 2 gives us 
additionally N 2 = 2d — 3 linearly independent combinations 
of Pm (the total number of second elements Hjfli in all rows), 
the third value 63 gives N3 = 2d — 5, and so on. 



TABLE I. The number of linearly independent POVM elements in- 
duced by m different quadrature measurements in a rf-dimensional 
Fock subspace. Bold font indicates IC POVMs. 
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An alternative way to get the same result is by noticing that 
after measuring the first quadrature, 2d — 1 linear combina- 
tions of pu are fixed and there remain d 2 — 1 — (2d — 1 ) = 
(d — 1 ) 2 — 1 free parameters entering the probability ©. This, 
effectively, lowers the dimension of the problem by one. Thus 
N 2 = 2(d-l)-l=2d-3,N 3 =2(d-2) -l = 2d-5, etc. 
The procedure can be repeated until the number of cuts m 
equals d, when Af reaches its maximum value of d 2 . One can 
draw an important conclusion from this discussion: by mea- 
suring each new quadrature, the prior information is updated 
in such a way that the dimension of the unknown part of the 
system becomes one unity lesser. 

In consequence, m different phase settings induce N = 
YL"k=\Nk linearly independent POVM elements; that is, 

!m(2d — m), m < d , 
(6) 
d 2 , m> d , 

which can be conveniently summarized as in Table J] This re- 
lation is the main achievement of this paper, as it provides a 
full description of the IC for the problem at hand. 

As we can see, to fully characterize a cZ-dimensional Fock 
subspace one needs d different quadratures. For a fixed m, the 
number N scales quadratically with d until d = m, due to the 
completeness of the generated POVMs. For d > m the size 
of POVM growth only linearly, hence in the asymptotic case 
of sufficiently large dimensionality the measurement becomes 
more and more incomplete. 

In this respect, we wish to highlight that our results do 
not rely on a particular reconstruction scheme, as in Ref. 10], 
but rather they emerge from the properties of the measure- 
ment itself and pervade all the continuous-variable detection 
schemes. 

Physical discussion. — The analysis presented so far has 
significant consequences, some of which will be explored 
here. Notice that the IC quadrature operators with d phases 
approximate the structure of mutually unbiased measure- 
ments. Without the Hilbert-space truncation, the projectors 
with the same phase setting are orthogonal, while for different 
phases they overlap. These complementary properties become 
modified in a finite-dimensional subspace, although they ap- 
proach their continuous-variable counterpart as long as the 
cutoff d is sufficiently large. 
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Loosely speaking, the orthogonal measurements bring al- 
ways a completely new bit of information, whereas the over- 
lapping ones verify the previously acquired information and 
add a new piece thereof. This is the essence of the tomogra- 
phy, where afterwards all the space is covered by an "oper- 
ator" mesh. Since the homodyne detection on finite dimen- 
sional subspaces is feasible, this gives an excellent opportu- 
nity to investigate the information gained in the process of the 
measurement. 

We emphasize that even a single quadrature actually gen- 
erates IC POVMs, provided that the signal space is properly 
matched to such a detection. For example, Table Q] shows 
that in a five-dimensional Fock subspace a single quadrature 
generates 9 independent POVM elements. These are IC in a 
suitably chosen three-dimensional subspace. Put differently, 
a single quadrature provides characterization of a qutrit in a 
five-dimensional space or two entangled qubits in an nine- 
dimensional space. Similarly, all Pauli matrices are induced 
by just two quadratures (say, x and p). Such kind of prior 
knowledge about the signal space can significantly reduce the 
number of measurements required for the complete state re- 
construction. A particularly natural situation is the case of 
symmetric/antisymmetric subspaces iflijl: a complete charac- 
terization of a two-qutrit state requires 9 different quadratures, 
whilst an antisymmetric two-qutrit state can be inferred with 
just one quadrature measurement. 

The previous results can be simplified in some cases where 
prior information is available. For instance, some elements of 
the density matrix can be known to vanish, and we can reach 
the complete POVM with less settings. Let us consider an 
analog of a squeezed vacuum of the form 

|¥) = c |0) +ci|4) + c 2 |8) + . . . + c k \4k) + .... (7) 

The coefficients q. gradually decrease with Ic, which makes 
possible to set them to zero at some threshold value = for 
k> d. The resulting state lives in a <i-dimensional subspace, 
so one may think that m = d phase settings would be needed 
to characterize that subspace. Surprisingly, the prior infor- 
mation about photons being generated in quadruples modi- 
fies the structure of the induced POVM and the number of 
linearly independent measurements is found to be different 
from that calculated before. Let us have a look at the case 
d = 3. The coefficients in (|4]i contain polynomial powers up 
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By writing explicitly the power expansion of the coefficients 
H H S = 256x 8 - 3584x 6 + 1344(k 4 - 13440x 2 + 1680 and 
H 4 H 4 = (16x 4 — 48x 2 + 12) 2 , we see that they are linearly in- 
dependent. Indeed, HoH% cannot be written as a linear com- 
bination of HqHq, HqH 4 and H 4 H 4 due to the lack of all 
needed polynomial powers, in contrast to the general case. 



Hence, the number of linearly independent POVM elements 
is much higher now: instead of N\ = 2d — 1, N2 = 2d — 3, 
N3 = 2d — 5, etc, we have the modified expressions N\ = 
Ni +N 3 +N 5 + ...,N%=N 2 ,NI=N 4 ,NZ=N 6 ,.... The num- 
ber of different phase settings for the scheme to be complete is 
m* = [d/2] + 1, where [] denotes the integer part. Hence, the 
family of states (O can be characterized with roughly half as 
many quadratures that would be required for a truncated Fock 
subspace of the same dimension. 

As our last contention, we insist in that proper matching 
conditions are essential l20ll . Particularly, without a deeper 
understanding, too dense binning of quadratures may be use- 
less, whereas quadrature phase may easily become undersam- 
pled. In other words, the number of quadrature bins multiplied 
by the number of phases settings does not, in general, equal 
the number of free parameters in the homodyne detection. 

Photon-number resolving detection. — To obtain our main 
result (O we have essentially relied on the polynomial rep- 
resentation of the quadrature ||3), namely that the probabil- 
ity distribution (01 can be written with help of polynomials 
up to a certain power, defined by dimensionality of the quan- 
tum state. This property holds also for many other situations; 
indeed, almost any physical measurement satisfies this prop- 
erty Il2tl . For example, an approximation of coherent states 

| a) = e~' a ! ' 2 £" = o a n /\fn\\n) by a finite sum, neglecting 
tiny contributions of sufficiently high n, is essentially the same 
situation as above: linearly independent POVM elements can 
be obtained from projections of a finite-dimensional state onto 
a properly matched basis. 

In the simplest scenario, the strategy can be photon- 
number-resolving detection. In fact, instead of we have 
now 

P(a,n) = |(„|a)| 2 = e -l«l 2 !^. (9) 

An obvious drawback of this choice is phase insensitivity, 
since coherent states with the same amplitude and different 
phases give identical statistics. To avoid this degeneracy, the 
photon-number detection can be preceded by a suitable dis- 
placement lfl9ll or even more complicated "steering" opera- 
tion. A detailed examination of this point is out of the scope 
of the present work. 

Concluding remarks. — We have addressed the fundamen- 
tal question of how many independent measurements can be 
generated by a given experimental setup. The analysis clari- 
fies the issue of necessary and sufficient complexity for tomo- 
graphical schemes. It turns out that the effective size of quan- 
tum tomography for any finite-dimensional system can be re- 
lated to the number of different quadratures observed, which 
illustrates an exciting link between discrete and continuous- 
variable systems. We stress that such a reconstruction cannot 
be accomplished without the appropriate delimitation of the 
reconstruction space, since the space touched by the observa- 
tion exceeds then the amount of available data. 

Our central result indicates that a single quadrature can 
be discretized up to 2d — 1 bins for a li-dimensional systen, 
This appears to be very similar to the classical Nyquist fre- 
quency in the Kotelnikov-Shannon coding theorem. By find- 
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ing a simple analytical rule we have provided a full charac- 
terization of informational completeness for a wide class of 
continuous-variable measurements and made a step towards a 
better understanding of quantum resources in quantum optics 
and quantum information processing. 

The analysis is applicable for arbitrary dimensionality and 
thus serves as a good alternative to other traditional ap- 
proaches based, e.g., on mutually unbiased bases and sym- 
metric IC POVMs. 

Finally, note that only data-acquisition issues have been ad- 
dressed here: the subsequent reconstruction must be capable 
to deal with generic nonequivalent detections. The mathe- 



matical method best suited for this purpose it the maximum 
likelihood estimation; however, even in this case, the contin- 
uous data are discretized due to the very nature of the mea- 
surement and the finite amount of memory and computational 
time available. We believe that our theory can be helpful in 
future optimizations of continuous-variable measurements. 

Financial support from the EU FP7 (Grant Q-ESSENCE), 
the Spanish DGI (Grant FIS201 1-26786), the UCM-BSCH 
program (Grant GR-920992), the Czech Technology Agency 
(Grant TEO 1020229), the Czech Ministry of Trade and Indus- 
try (Grant FR-TI 1/384) and the IGA of Palacky University 
(Grant PRF- 20 12-005) is acknowledged. 



[1] E. Prugovecki, Int. J. Theor. Phys. 16, 321 (1977); P. Busch and 

P. J. Lahti, Found. Phys. 19, 633 (1989). 
[2] M. G. A. Paris and J. Rehacek, eds., Quantum State Estimation, 

Lect. Not. Phys., Vol. 649 (Springer, Berlin, 2004). 
[3] G. M. D. Ariano, P. Perinotti, and M. F. Sacchi, J. Opt. B 6, 

S487 (2004); S. T. Flammia, A. Silberfarb, and C. M. Caves, 

Found. Phys. 35, 1985 (2005); S. Weigert, Int. J. Mod. Phys. B 

20, 1942 (2006); T. Durt, Open Sys. Inf. Dyn. 13, 403 (2006). 
[4] K. Vogel and H. Risken, Phys. Rev. A 40, 2847 (1989). 
[5] D. T. Smithey, M. Beck, M. G. Raymer, and A. Faridani, Phys. 

Rev. Lett. 70, 1244(1993). 
[6] J. Rehacek, S. Olivares, D. Mogilevtsev, Z. Hradil, M. G. A. 

Paris, S. Fornaro, V. D'Auria, A. Porzio, and S. Solimeno, 

Phys. Rev. A 79, 032111 (2009). 
[7] U. Leonhardt and M. Munroe, Phys. Rev. A 54, 3682 (1996); 

U. Leonhardt, J. Mod. Opt. 44, 2271 (1997). 
[8] A. V. Oppenheim, A. S. Willsky, and S. Hamid, Signals and 

Systems, 2nd ed. (Prentice Hall, New lersey, 2010). 
[9] J. Schwinger, Proc. Natl. Acad. Sci. USA 46, 570 (1960); W. K. 

Wootters and B. D. Fields, Ann. Phys. 191, 363 (1989); A. B. 

Klimov, J. L. Romero, G. Bjork, and L. L. Sanchez-Soto, 

J. Phys. A 40, 3987 (2007); A. B. Klimov, D. Sych, L. L. 

Sanchez-Soto, and G. Leuchs, Phys. Rev. A 79, 052101 (2009); 

T. Durt, B.-G. Englert, I. Bengtsson, and K. Zyczkowski, Int. 

J. Quantum Inf. 8, 533 (2010). 



[10] R. L. Stratonovich, J. Stochastics 1, 87 (1973). 

[11] B. A. Grishanin, Tekn. Kibernetika 11, 127 (1973). 

[12] C. W. Helstrom, Quantum Detection and Estimation Theory, 
(Academic, New York, 1976). 

[13] A. S. Holevo, Probabilistic and Statistical Aspects of Quantum 
Theory, 2nd ed. (North Holland, Amsterdam, 2003). 

[14] V. B. Braginski and F. Y. Khalili, Quantum Measurement (Ox- 
ford University Press, Oxford, 1992). 

[15] C. W. Helstrom, J. W. S. Liu, and I. P. Gordon, Proc IEEE 58, 
1578 (1970). 

[16] J. M. Renes, R. Blume-Kohout, A. I. Scott, and C. M. Caves, 
I. Math. Phys. 45, 2171 (2004); D. M. Appleby, S. T. Flam- 
mia, and C. A. Fuchs, ibid. 52, 022202 (2011); R. Salazar, 
D. Goyeneche, A. Delgado, and C. Saavedra, Phys. Lett. A 
376, 325 (2012). 

[17] J. Rehacek, Z. Hradil, Z. Bouchal, R. Celechovsky, I. Rigas, 
and L. L. Sanchez-Soto, Phys. Rev. Lett. 103, 250402 (2009). 

[18] D. Sych and G. Leuchs, New I. Phys. 11, 013006 (2009); 11, 
013006 (2009). 

[19] C. Wittmann, U. L. Andersen, M. Takeoka, D. Sych, and 
G. Leuchs, Phys. Rev. Lett. 104, 100505 (2010); Phys. Rev. 
A 81, 062338 (2010). 

[20] A. Ourjoumtsev, R. Tualle-Brouri, I. Laurat, and P. Grangier, 
Science 312, 83 (2006). 



